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Abstract 



> 

in 

\f~} ' We compute the Taylor expansion for the luminosity distance for an observer at the center of a 

spherically symmetric matter inhomogeneity with a non vanishing cosmological constant. Defining 

(N ■ 

apparent the value of the cosmological constant estimated by an observer who observes the same 
luminosity distance but assumes homogeneity in its estimation of the cosmological parameters, we 
derive a relation between the apparent and the true value, i.e. the one obtained taking into account 
the inhomogeneity. We can consider the assumption of being at the center of a spherically symmetric 
inhomogeneous matter distribution equivalent to calculating the monopole contribution of the large 
scale structure around us, which we expect to be the dominant one, because of other observations 
such as the cosmic microwave background radiation supporting a high level of isotropy. 
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I. INTRODUCTION 

Luminosity distance measurements [l -f| and the WMAP measurements 7, 8] of cosmic mi- 
crowave background (CMB) interpreted in the context of standard FLRW cosmological models 
strongly support a dominant dark energy component. 

As an alternative to dark energy, some of attenton has been devoted to the possib ility o: 
interpreting the same observational data with inhomogeneous comsological model 9l.ll0l.ll2l-u 

Our approach will be different, since we will consider a Universe with a cosmological constant 
and some local large scale inhomogeneity modeled by a ALTB solution 30]. For simplicity we 
will also assume that we are located at its center. In this regard this can be considered a 
first attempt to model local large scale inhomogeneities in the presence of the cosmological 
constant or, more in general, dark energy. Given the spherical symmetry of the LTB solution 
and the assumption to be located at the center our calculation can be interpreted as the 
monopole contribution of the large inhomogeneities which surround us. Since we know from 
other observations such as CMB radiation that the Universe appears to be highly isotropic, we 
can safely assume that the monopole contribution we calculate should also be the dominant 
one, making our results even more relevant. We fist calculate the null radial geodesies for a 
central observer and then use it to obtain the luminosity distance in a LTB model with non 
vanishing comsological constant. We then compare it to that of ACDM models, finding the 
relation between the two different cosmological constants appearing in the two models, where 
we call apparent the one in the ACDM and true the one in ALTB. 

As an applicaton we then calculate the local Taylor expansion of the function defining an 
LTB model which is able to accound for a 5% difference between the apparent and the true 
value of the cosmological constant. 



II. LTB SOLUTION WITH A COSMOLOGICAL CONSTANT 



The LTB solution can be written as 



31 



33( as 



ds 2 



-dt 2 + 



(R, r f dr 2 



+ R 2 dQ 2 



1) 



l + 2£(r 

where R is a function of the time coordinate t and the radial coordinate r, E(r) is an arbitrary 
function of r, and R r = d r R(t,r). The Einstein equations with dust and a cosmological 



2 



constant give 

(R\ 2 2E(r) 2M(r) A 

~ f? 2 + i?3 + 3' lJ 

^.O = (3) 

with M(r) being an arbitrary function of r, i? = d t R(t, r) and c = 87rG = 1 is assumed 
throughout the paper. 

Eq.(|2J) is obtained by integration respect to time of one of the Einstein's field equations, 
and the function M(r) plays in fact the role of constant of interation. Since Eq. fl2]) contains 
partial derivatives respect to time only, its general solution can be obtained from the FLRW 
equivalent solution by making every constant in the latter one an arbitrary function of r. 

The general analytical solution for a FLRW model with dust and cosmological constant was 



obtained by Edwards |34| in terms of elliptic functions. By an appropriate choice of variables 
and coordinates, we may extend it to the LTB case thanks to the spherical symmetry of both 
LTB and FLRW models, and to the fact that dust follows geodesies without being affected 
by adjacent regions. An anaytical solution can be found by introducing a new coordinate 
7] = rj(t, r) and a variable a by 



' drj\ r _ 1 
~dt) ~R = a 



(4) 



and new functions by 



Then Eq. (T5]) becomes 



_ 6M(r) _ 2E(r) . 
P^ r ) = ' fc ( r ) = -T- ■ ( 5 ) 



(da\ 2 2 p (r) A 4 

= - Hr)a + — a+ 3 a ' (6) 

where a is now regarded as a function of rj and r, a = a(rj,r). It should be noted that the 

coordinate rj, which is a generalization of the conformal time in a homogeneous FLRW universe, 

has been only implicitly defined by Eq. Q. The actual relation between t and r\ can be obtained 

by integration once a(r], r) is known: 

rv 

t(r),r)= / a(x,r)dx + t b (r) , (7) 
Jo 



35|]. 



which can be computed analytically, and involve elliptic integrals of the third kind 

The function t_e(r) plays the role of constant of integration, and is an arbitrary function of r, 
sometime called bang function, since by construction at time t = we have a(i&(r),r) = 0, 



and correspond to the fact that the big bang initial singularity can happen at different times 
at different positions from the center in a LTB space. In the rest of this paper we will assume 
homogeneous bang, i.e. we will set 

t b (r) = 0. (8) 
Inspired by the construction of the solution for the FLRW case get: 

<V,r) = — j- 7TT , (9) 

where (f)(x; g 2 , g 3 ) is the Weierstrass elliptic function satisfying the differential equation 



2 



dx 



920 - 9s , (10) 



and 



4 4 
a = Po{r) , 92 = ^k(r) 2 , g 3 = — (2A;(r) 3 - Ap (r) 2 ) . 



'ir, 



3 v , , 27 

In this paper we will choose the so called FLRW gauge, i.e. the coordinate system in which 
Po(r) is constant. 



It is convenient to re-write the solution in terms of dimensionless quantities 



33| : 



k(r) 
V 
Po 
A 

a(rj,r) 



(a H ) 2 K(r) , 
T{a H )-\ 

oqO n 3 tt2 

a{T{a H )~ 1 , r) = a(T,r) , 



(12) 
(13) 
(14) 
(15) 
(16) 
(17) 



to obtain 



a(T, r) 

92(r) 
9s(r) 



M 



K(r) + 12<P(T,g 2 (r),g 3 (r)) 
K(r) 2 



12 

^- (2K(r f - 270 A (0° M ; 



(19) 
(20) 
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or after multiplying every term by (a^Ho) 2 and using the original dimensionful quantities 



a(V,r) k(r) + 12^, g 2 (r),g 3 (r)) fi(T ' r) ' 
0(77,r) = (f)(r](a H ),r)(a H ) 2 = (p(T,r)(a H ) 2 . 



(21) 
(22) 



In this form is an arbitrary scale which we can set equal to the observed value and to the 
Hq TB by appropriately setting the value of T as explained in more details in [371 ] . 



III. GEODESIC EQUATIONS AND LUMINOSITY DISTANCE 



We use the same method developed in 



361 ] to solve the null geodesic equation written in 
terms of the coordinates (r),r). Rather than integrating differential equations numerically, we 
perform a local expansion of the solution around z = corresponding to the point (t ,0), or 
equivalently (770,0), where to = t(rjo,0). 

The luminosity distance for a central observer in the LTB space-time as a function of the 
redshift z is 

D L (z) = (l + z) 2 R(t(z),r(z)) = (l + z) 2 r(z)a( V (z),r(z)) , (23) 

where (t(z),r(zfj or ((r](z),r(z)) is the solution of the radial geodesic equation as a function 
of z. The equation for radial null geodesies can be easily obtained in the coordinates (t, r) 

dt R,r{t,r) 



dr 



1 + 2E(r) 



where t = T(r) is the time coordinate along the nul 
coordinate r. In terms of z, Eq. (1241) takes the form 38 



(24) 



radial geodesic function of the 



dr 
dz 
dt 
dz 



l + 2E(r(z)) 



[l + z)R r [r(z),t(z)} 
R r [r(z),t(z)} 



(l + z)R, r [r(z),t(z)} 
These equations are derived from the definition of redshift and by 
short time interval along the null geodesic T(r). It can be shown 
(77, r) eqs. (|25|) take the form: 



(25) 

owing the evolution of a 



361 ] that in the coordinates 



ol 



F(r),r) 



d r (a(r), r)r) — a d v (a(r), r)r) d r t(r], r) 



(28) 



where 

R r 1 
T + 2E(r) yjl - k(r)r 2 

It is important to observe that the functions p, q, F have explicit analytical forms, making it 
particularly useful do derive anaytical results. 



IV. FORMULA FOR THE LUMINOSITY DISTANCE 

In order to obtain the redshift expansion of the luminosity distance we need to use the 
following: 

k(r) = (a H ) 2 K(r) = K + K x r + K 2 r 2 + .. (29) 

t(r),r) = b {r}) + h^r + b 2 (rj)r 2 + ■■ (30) 

We will use eq.(7) to obtain the above expansion for t(rj,r) from the exact solution for 



a(rj,r). Following the same approach given in [30] , we can find a local Taylor expansion in 
red-shift for the geodesies equations, and then calculate the luminosity distance. The general 
expression is rather cumbersome, so will report here only the result assuming K = 0, which 
is still showing the general nature of the effect. We will expand the solution of the geodesic 
equations according to: 

r(z) = riz + r 2 z 2 + ... (31) 
T)(z) = 7]\Z + r] 2 z 2 + ... (32) 
K(z) = K lZ + K 2 z 2 + ... (33) 

(34) 

After substituting the above expansion in the geodesies equation we can reduce the solution of 
the differential equations to the solution of a system of algebraic equations for the expansion 
coefficient. 

For the geodesies equations we get: 

iTi(To-l)To + 3n M 



3a H Q M 
1 



m 36a H n A n 2 



3n A n M (9n 2 M - ak 2 (-i + t )t ) + 3x x q a o m (-4 + (4 - 9O M )r + 

Mi 

+(-4 + 9Q M )r 2 ) + ^(2^(2 + (-4 + 3Q M )T - 6(-l + n M )T 2 + 3(-l + Q M )Tl 



-4(-l + T )WZ + fi M (-4 + T + 3T WZ)) 
1 



r 3 



a H ' 



12a H Q M 

'2 



(9Q 2 M + (-4 + (4 - 6Q M )T + (-4 + 6ft M )T 2 ) 



72a H n A ni 



M 



2^(6(3^ - 4Q M + 1)T 4 - 12(3^ - 4Q M + 1)T 3 + 3(6Q 2 M - 13ft 



Ai 



(36) 
(37) 
(38) 
4)T 2 + 



+2(9ft M - 4)T + 4) + 3fti f T (3T Co + T - 4) - 2ft M (9T 2 Co + T 2 - 6T Co - 4T 
6Co - 2) + 8(T 2 - T + l)Co) + 12/^ft A ft 2 f ((9ftA/ - 8)T 2 + (8 - 9tt M )T - 5) + 
+3n A n M (K 2 ((8 - 12ft M )T 2 + 4(3ft A/ - 2)T + 8) + 3(9ft M - 4)ft 2 f )) 

After substituting in the formula for the luminosity distance and expanding we finally get: 



(39) 



D 



ALTS, 
L 1 



D 
D 



ALTB 



ALTB 
2 

ALTB 



(1 + z) 2 r(z)a 
1 

1 



ALTB, 



( V (z),r(z)) = Dt L1B z + D^^z 2 + D 



\ALTB 3 



4# r 



(-4 + 3ft M + 2K 1 (-l + ro)r ) 



ir 2 (2ft A T ((6ft M - 5)T 3 - 2(6ft M - 5)T 2 + 6(ft M 



+n M T (3T ( + T - 4) - 4(T 2 Co - T ( + Co - 1)) 
+4ir 1 ft A ft M ((9ft M - 8)T 2 + (8 - 9n M )T - 2) + 
+3ft A ftA/(-4K 2 (T " l)^o + 9J& - l0tt M ] 
where we used the Einstein equation at the center (77 = rj , r = 0) 

1 = ftfc(0) + Sljw + ^A = -ifo + ftjVf + ^A, 
Mr") 

^A 



ft fe (0) + Q M + tt A 
k(r) 

Po 
A 



and To = ?yo(ao-f/o) is determined numerically by imposing the conditions 37] 



a(£, r 
0(770,0) = a , 



d t a(t,r) d v a(r],r) a'(T,r) 
- {a H )- 



0(77, r) 



a(T, r) 5 



# LTB (%,0) = flo- 
In the above formulae we have also defined 

Co = Cfab,02(o),0 3 (o)) 



(40) 
(41) 



(42) 
1)T + 2) + 



(43) 

(44) 
(45) 

(46) 
(47) 

(48) 
(49) 
(50) 

(51) 



where ( is the Weierstrass Zeta Function satisfying the equation 

d((z,92(r),g 3 (r)) 

^ = -nz,92{r),g 3 {r)) , (52) 

The presence of this last function in the formulae obtained above is due to the fact that the 
function t(j],r) which enters the geodesies equation is the integral of a(r],r), and since this 
depends on <f)(z), its integral will depend on ((z). In the case of a LTB solution without 
a cosmological constant this integral can be performed without the introduction of a new 
function, while in this case it requires the introduction of Co i n the final formula. 

The procedure to reduce the analytical formula to this form is rather complicated since it 
involves to express wherever possible all the intermediate expressions in terms of physically 
meaningful quantities and to use the properties of the elliptic functions. We give more details 
about it in the appendix. 
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FIG. 1: The percentual error A = 100 Tav ^°i TB Num between the numerically computed DjfiJ^ and 

Taylor 

the Taylor third order expansion D^^ lor is plotted as a function of the redshift is plotted for the 
LTB solution corresponding to K(r) = K\r + K2r 2 . The coefficients K\, K2 are obtained by matching 
the third order expansion of the luminosity distance of a ACDM model using the Taylor expansion 
obtained in this paper, but with a smaller value of the cosmological constant &j( rB = 95%flj^ DM . 
As it can be seen the Taylor aproximation is quite accurate at low red-shift. 
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V. CALCULATING D L (z) FOR KCDM MODELS. 



The metric of a ACDM model is the FLRW metric, a special case of LTB solution, where : 

Po(r) oc const, (53) 

k{r) = 0, (54) 

t b (r) = 0, (55) 

a(t,r) = a(t). (56) 

We will calculate independently the expansion of the luminosity distance and the redshift 
spherical shell mass for the case of a flat AC DM, to clearly show the meaning of our notation, 
and in particular the distinction between f2^ pp and Vt\ ue . We can also use these formulas to 
check the results derived before, since in absence of inhomogeneities they should coincide. 
One of the Einstein equation can be expressed as: 



H ACDM (z) = # y (1 - fiD (7) 3 + n T = #oV(l - P )(l + zf + «T- (57) 
We can then calculate the luminosity distance using the following relation, which is only valid 
assuming flatness: 

Df DM (z) = (l + z) [ H J D Z ' M{zl) = Df DM z + D^™ Z 2 + d acvm z s + m 
From which we can get: 

Dt CDM = ^ , (59) 

jjkCDM = ^+1 (6Q) 

We can check the consistency between these formulae and the ones derived in the case of LTB 
by setting: 

Kt = K = , (61) 

(62) 

which corresponds to the case in which fi A pp = f2^ ue . 

VI. CORRECTION TO THE APPARENT VALUE OF THE COSMOLOGICAL CON- 
STANT 

As explained in [40], the matching of the first two order terms of the Taylor expansion of 
the luminosity distance of a ALTB and a AC DM model are sufficient to obtain the leading 



order relation between the apparent and true value of the cosmo logical constant. The matching 
of higher orders terms gives in fact relations necessary to set the other coefficients of the 
inhomogeneity profile, but do not affect the relation among Q app and Q true . After calculating 
the first two terms of the redshift expansion of the luminosity distance for ALTB and AC DM 
model we can now look for the ALTB models which give the same theoretical prediction the 
best fit ACDM model. 

From the above relations we can derive : 



K, 



12(0^ - l)Qi(T - 1] 
+Co - 1) - 4T (2J4™ e + Co - 1) + 4(Co - 1)) + (64) 

-4(ftf e - l)^™ e (3T 2 - 3T + 2)) . (65) 
As expected the relation above reduces to 

in the limit in which there is no inhomogeneity, i.e. when K\ = 0. We should remember that 
all the formulae above are derived under the assumption that K = 0, and that the value 
of Tq depends on Q 1 ™ 6 and K through the procedure which was explained in the previous 
sections. A more general formula can be derived also the case in which Kq is not zero, but 
given ts rather complicated form we do not report it here, since the one derived above already 
captures the essence of the effect of the large scale inhomogeneity on the apparent value of 
the cosmological constant. We can see in fact that the correction can be positive or negative 
depending on the sign of K\, which is in qualitative agreement with previous studies of the 
effects of inhomogeneities on the effective equations of state of a ALTB cosmological model 
H. 

As an application we assume a 5% difference between the apparent and the true value of the 
cosmological constant and compute Ki, Ki from the matching of the coefficients of the Taylor 
expansion of the luminosity distance, obtaining: 

K x = -0.007462, (67) 
K 2 = 0.007037. (68) 

This gives the ALTB model able to mimic locally such a difference in the value of Vt l ^ ue and Vt a £ p . 
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The accuracy of the formula is then compared to the numerical calculation of the luminosity 
distance, showing a good agreement. 

VII. CONCLUSIONS 

We have derived the correction due to large scale structure to the value of the apparent 
cosmological constant inferred from low redshift supernovae observations. This analytical cal- 
culation shows how the presence of a local inhomogeneity can affect the estimation of the value 
of cosmological parameters, such as Q\. This effects should be properly taken into account 
both theoretically and observationally. It is important to underline here that we do not need a 
large void as normally assumed in previous studies of LTB models in a cosmological context. 
Even a small inhomogeneity could in fact be important. The main advantage of the analytical 
approach we adopted is that it allows to understand the effects in a model independent way, 
without making any assumption about the nature of the inhomogeneity, which could be either a 
local under or over density. The formula we derived predicts in fact that the cosmological could 
be both bigger or smaller than the apparent value obtained by assuming large scale homogene- 
ity. In the future it will also be interesting to extend the same analysis to other observables 
such as barionic acoustic oscillations (BAO) or the cosmic microwave background radiation 
(CMBR), and we will report about this in separate papers. Another direction in which the 
present work could be extended is modeling the local inhomogeneity in a more general way, for 
example considering not spherically symmetric solutions. From this point of view our calcula- 
tion could be considered the monopole contribution to the general effect due to a local large 
scale inhomogeneity of arbitrary shape. Given the high level of isotropy of the Universe shown 
by other observations such as the CMB radiation, we can expect nevertheless the monopole 
contribution we calculated to be the dominant one. 

Appendix A: derivation of the analytical formulae 

In order to obtain the formula for Dl(z) in the form in which we reported it in the previous 
sections, we need to apply several simplifying procedures. From the definition of a and H we 
can get: 

0o = #7o,<? 2 (O),S 3 (O)) = P ° ~ a ° k0 , (1) 

Sa 
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M7,<72(0),g 3 (0)) 
dr] 



2H oPo 
3 



(2) 



We can then substitute the above expressions everywhere {4>o,4>b} appear, which is the reason 
why they are not present in the formulae obtained in the previous sections. 

Another useful relation is the one which can be obtained from the differential equation 
defining the Weierstrass elliptic function eq. ( TTOj) : 



which can be used for example to obtain T from eq. f)49ll50p . since it allows to determine 0o 
analytically and then Tq numerically by solving the equation : 
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